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Iron-based superconductivity develops near an antiferromagnetic order and out of a bad metal 
normal state, which has been interpreted as originating from a proximate Mott transition. Whether 
an actual Mott insulator can be realized in the phase diagram of the iron pnictides remains an open 
question. Here we use transport, transmission electron microscopy. X-ray absorption spectroscopy, 
and neutron scattering to demonstrate that NaFei_a,Cua;As near x ~ 0.5 exhibits real space Fe and 
Cu ordering, and are antiferromagnetic insulators with the insulating behavior persisting above the 
Neel temperature, indicative of a Mott insulator. Upon decreasing x from 0.5, the antiferromagnetic 
ordered moment continuously decreases, yielding to superconductivity around x = 0.05. Our discovery 
of a Mott insulating state in NaFei_a,Cu,5As thus makes it the only known Fe-based material in which 
superconductivity can be smoothly connected to the Mott insulating state, highlighting the important 
role of electron correlations in the high-Tc superconductivity. 


INTRODUCTION 


At the heart of understanding the physics of the iron-based superconductors is the interplay of superconductivity, 
magnetism and bad-metal behavior [iHg , a key question is whether superconductivity emerges due to strong electronic 
correlations iSIl or nested Fermi surfaces [l2l - [ll| . Superconductivity occurs in the vicinity of antiferromagnetic (AF) 
order, both in the iron pnictides and iron chalcogenides [m. In addition, the normal state has a very large room- 
temperature resistivity, which reaches the loffe-Mott-Regel limit [las- This bad-metal behavior has been attributed 
to the proximity to a Mott transition [M, with the Coulomb repulsion of the multi-orbital 3d electrons of the Fe 
ions being close to the threshold for electronic localization. In the iron chalcogenide family, several compounds have 
been found to be AF and insulating with characteristic features of a Mott insulator lB“18|. There is also evidence 
for an orbital-selective Mott phase in Aa;Fe 2 -ySe 2 {A = K, Rb) [l^. However, in these iron chalcogenide materials, 
one cannot continuously tune the AF Mott insulating state into a superconductor. On the other hand, the iron 
pnictide NaFei_a;Cua;As is a plausible candidate for a Mott insulator based on transport and scanning tunneling 
microscopy (STM) measurements near x = 0.3 (2^, [2lj, however the insulating behavior may also be induced by 
Anderson localization via Cu-disorder. 

The parent compounds of iron pnictide superconductors such as AFe 2 As 2 (A =Ba, Sr, Ca) and NaFeAs have crystal 
structures shown in Fig. 1(a) and 1(b), respectively 0. They exhibit a tetragonal-to-orthorhombic structural phase 
transition at temperature Tg, followed by a paramagnetic to AF phase transition at Tn (Ts > 7n) with a collinear 
magnetic structure where the spins are aligned antiferromagnetically along the a-axis of the orthorhombic lattice [Fig. 
1(c)] [ 2 ^. When Fe in AFe 2 As 2 is replaced by Cu to form ACU 2 AS 2 241, the Cu atoms have a nonmagnetic 3(i^® 
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electronic configuration with Cu^'*' and a filled d shell due to the presence of a covalent bond between the As atoms 
within the same unit cell [As]“^ = [As-As]“^/2 shown in the shaded As-As bond of Fig. 1(a) (25l - E8| . as predicted 
by band structure calculations [ 2 ^. Since the crystal structure of NaFei_ 2 ;Cua;As does not allow a similar covalent 
bond [Fig. 1(b)] [^, it would be interesting to explore the magnetic state of NaFei-a^Cu^^As in the heavily Cu-doped 
regime. From transport measurements on single crystals of NaFei_a;Cua:As with x < 0.3, it was found that the in¬ 
plane resistivity of the system becomes insulating-like for x > 0.11 [ 2 ^. STM showed that local density-of-states for 
NaFei_a,Cua;As with x = 0.3 resembles electron-doped Mott insulators [2l|, although it is unclear what the undoped 
Mott insulating state is or whether it exhibits AF order. 

Using transport, neutron scattering, transmission electron microscopy (TEM), X-ray absorption spectroscopy 
(XAS), and resonant inelastic X-ray scattering (RIXS), we demonstrate that heavily Cu-doped NaFei_a:Cua;As ex¬ 
hibits Fe and Cu ordering and becomes an antiferromagnetic insulator when x approaches 0.5. The insulating behavior 
persists into the paramagnetic state, providing strong evidence for a Mott insulating state. This conclusion is cor¬ 
roborated by our theoretical calculations based on the 17(1) slave-spin approach [s^. Our calculations demonstrate 
enhanced correlations and a Mott localization from the combineed effect of a bandwidth reduction, which results from 
a Cu-site blockage of the kinetic motion of the Fe 3d electrons, and a hole doping from Sd® to 3d®, both of which 
are made possible by the fully occupied Cu 3d shell BiSil l30l |. Upon decreasing x from x ~ 0.5, the correlation 
lengths of Fe and Cu ordering and magnetic ordering, as well as the ordered magnetic moment continuously decrease, 
connecting smoothly with the superconducting phase in NaFei_xCua;As appearing near x = 0.05, highlighting the 
role of electronic correlations in iron pnictides. 


RESULTS 


Resistivity measurements 

Single crystals of NaFei-xCu^As were prepared using the self-flux method. Transport measurements were car¬ 
ried out using a commercial physical property measurement system (PPMS) with the standard four-probe method 
[Methods section]. Figure 1(e) shows temperature dependence of the in-plane resistivity for the x = 0.016 sample. 
In addition to superconductivit y a t Tc = 11 K, the normal state resistivity p is ^ 0.4 mil • cm at room temperature, 
consistent with previous work [20|. Figure 1(f) plots the temperature dependence of the resistivity on a log scale 
for samples with x = 0.18, 0.38, 0.44 and 0.48. Resistivity in the x = 0.18 sample exhibits insulating-like behavior 
consistent with earlier work [20|. For the x > 0.39 samples, resistivity further increases by an order of magnitude 
compared to the x = 0.18 sample, signifying that electrons become much more localized in NaFei_xCua;As when x 
approaches 0.5. 


Neutron scattering results on NaFeo.s6Cuo.44As 


Unpolarized and polarized elastic neutron scattering were carried out on NaFei_ 2 ,Cu 2 ;As single crystal samples 
[Methods section]. For these measurements, we use the orthorhombic unit cell notation suitable for the AF ordered 
state of NaFeAs [13, 23], and define momentum transfer Q in three-dimensional reciprocal space in as Q = 
Hsl* -|- Kh* 4- Ac*, where H, K, and L are Miller indices and a* = a 27 r/a, b* = b 27 r/ 6 , c* = c2TTjc [Fig. 1(c)]. 
Single crystals are aligned in the {H,0,L) and {H,H,L) scattering planes in these measurements. In the (AT, 0,A) 
scattering geometry, the collinear magnetic structure in NaFeAs [Fi g. l (c) | gives magnetic Bragg peaks below Tn at 
Qaf = [H, 0, A), where H = 1, 3, • • • and A = 0.5,1.5, • • •, positions [221 l23l[. 

Figure 2(c) shows typical elastic scans for NaFeo. 56 Cuo. 44 As along the [iA, 0,0.5] direction at T = 3.6 K and 300 
K . Because of twinning in the orthorhombic state, this is equivalent to elastic scans along the [0,77,0.5] direction. 
While the scattering has a clear peak centered around Q = (1, 0,0.5)/(0,1, 0.5) at room temperature, the scattering 
is enhanced dramatically on cooling to 3.6 K, suggesting the presence of static AF order. The temperature difference 
plot between 3.6 K and 300 K in Fig. 2(d) indicates that the low-temperature intensity gain is essentially instrumental 
resolution limited (horizontal bar). Figure 2(e) shows the temperature difference along the [1, 0, A]/[0,1, A] direction 
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(a) BaCUjASj (b) NaFe, ^Cu As (c) NaFeAs (d) NaFe^ jCu^ jAs 





Temperature (K) Temperature (K) 

Fe 3d5 Fe 3c^ 


Mott insulator Bad metal 



FIG. 1: Summary of transport and neutron scattering results, (a) The crystal structure of j 4 Cu 2 As 2 in the tetragonal 
state, where ax is the in-plane lattice parameter. The As-As covalent bonding distance is ~2.5 A. (b) The crystal structure of 
NaFei_a;Cua;As, where similar As-As covalent bonding is not possible, (c) The collinear magnetic structure of NaFeAs, where 
the AF order and moment direction are along the orthorhombic a axis . Only Fe atoms are plotted in the figure for 

clarity, (d) Real space structure and spin arrangements of NaFeo. 5 Cuo. 5 As in the AF orthorhombic unit cell similar to NaFeAs. 
In NaFeo. 5 Cuo. 5 As, since Fe and Cu form stripes, Na and As also shift slightly from their high symmetry positions, (e) In-plane 
resistivity for x = 0.016 sample where bulk superconductivity occurs below Tc = 11 K [ 2 ^. (f) Temperature dependence of the 
in-plane resistivity for x = 0.18, 0.39, 0.44 and 0.48 samples. The vertical arrow marks the position of Tn for the x = 0.44 sample 
determined from neutron scattering, (g) Evolution of ordered moment with doping in NaFei-xCux As. The two ordered phases 
are separated by the superconducting dome marked as SC, with no magnetic order for x = 0.016 near optimal superconductivity. 
The right and left insets show the in-plane magnetic structures of NaFeAs and the new AF insulating phase, respectively. For 
X < 0.05, the phase diagram from Ref. [ 2 ^ is plotted. Vertical error bars are from least-square fits (1 s. d.). 


which apparently is not resolution limited. The temperature dependence of the scattering at (1,0,0.5)/(0,1, 0.5) 
plotted in Fig. 2(f) reveals clear evidence of magnetic ordering below Tn ~ 200 K, suggesting that the small peak 
observed at room temperature in Fig. 2(c) occurs in the paramagnetic phase and is thus of structural (super-lattice) 
origin induced by Cu substitution, since such a peak is forbidden for NaFeAs. 


Fe and Cu ordering in NaFei-^Cua, with x « 0.5 

To conclusively determine the origin of super-lattice peaks and the crystal structure of NaFeo. 56 Cuo. 44 As, we have 
carried out high-resolution TEM measurements. The inset in Fig. 2(a) shows an image of the sample, and the electron 
diffraction patterns were collected from areas about 100 nm in diameter within a single domain of the sample [see the 
arrow in the inset of Fig. 2(a)]. A typical diffraction pattern along the [001] zone axis at room temperature is shown 
in Fig. 2(a). While we see clear super-lattice reflections at H = 1,3, ■ ■ ■ positions along the [H, 0,0] direction, they are 
absent at the K = 1,3, • • • positions along the [0, K, 0] direction [Fig. 2(b)]. This means that the crystal structure of 
NaFeo. 56 Cuo, 44 As is orthorhombic and obeys two-fold rotational symmetry. From this information and from intensities 
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FIG. 2: TEM and neutron scattering results on the structure and magnetic order for NaFeo.56Cuo.44As. (a) 

Using 200 keV incident electrons, an electron diffraction pattern in the [H,K] plane was obtained from a NaFeo. 56 Cuo. 44 As 
particle with its TEM real-space image shown in the inset. The scale bar is 1/mi. (b) Cuts of (a) along the [H, 0] and [0, K] 
directions, (c) Unpolarized neutron diffraction scans along the [f/, 0, 0.5]/[0, A, 0.5] direction for the x = 0.44 sample at 300 
K and 3.6 K. The peak at 300 K is a super-lattice peak arising from Fe-Cu ordering. Temperature difference plots between 
3.6 K and 300 K along the (d) [A, 0, 0.5]/[0, A, 0.5] and (e) [l,0,L]/[0, 1,1/] directions. Solid lines are Gaussian fits and the 
horizontal bars indicate instrumental resolution, (f) Temperature dependence of the scattering at Q = (1,0, 0.5)/(0,1,0.5) 
shows Tn « 200 K. (g) Neutron polarization analysis of the magnetic Bragg peak at (1, 0, 0.5)/(0, 1, 0.5). Neutron SF and 
NSF cross sections are measured at 2 K and 240 K. The peak at 240 K is nonmagnetic nuclear scattering giving rise to NSF 
scattering, (h) Similar scans for the (1,1, 0) peak. The inset shows positions of these two peaks in reciprocal space. All vertical 
error bars represent statistical error (1 s. d.). 


of the super-lattice peaks in {H, 0, L) scattering plane, we conclude that Fe and Cu atoms in NaFeo. 5 Cuo. 5 As (which is 
approximated by NaFeo. 56 Cuo. 44 As) form a real space stripe-like ordered structure in Fig. 1(d) in the space group Ibam 
[Supplementary Note 1 and Supplementary Table 1]. This conclusion is further corroborated by single crystal neutron 
diffraction measurements [Supplementary Table 2]. The Fe-Cu ordering is a structural analogue of the magnetic order 
in NaFeAs and structural super-lattice peaks occur at the same positions as magnetic peaks in NaFeAs. Fe and Cu 
ordering in NaFei-^CuajAs affects the intensities of Bragg peaks already present in NaFeAs very little [Supplementary 
Note 1 and Supplementary Figure 1], and the induced super-lattice are relatively weak. Aside from ordering of Fe and 
Cu, other aspects of the structure of NaFei-^Cu^jAs remain identical to that of NaFeAs, corroborated by the similar 
neutron powder diffraction patterns [Supplementary Note 1, Supplementary Figure 2 and Supplementary Table 3]. 
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FIG. 3: Cu-doping evolution of the magnetic order in NaFei_a;Cua;As. (a) Comparison of wave vector scans along 
the [FI, 0,0.5]/[0, A, 0.5] direction for NaFei-a^Cua, As single crystals with x = 0.18,0.39, and 0.44. The data are normalized 
to the (2,0,0) nuclear Bragg peak. Note that intensity in the x = 0.18 and 0.39 samples are multiplied by 200 and 5 times, 
respectively, (b) Similar scans along the [1, 0, F/]/[0,1, L\ direction, (c) Cu-doping evolution of the spin-spin correlation length. 
The correlation length along (1,0,0) for the x = 0.44 sample is resolution limited, (d) Cu-doping evolution of the magnetic 
order parameter for NaFei-xCu^As. All error bars represent statistical error (1 s. d.). 


Magnetic structure of NaFei-ajCua,As with x « 0.5 


Having established the real space Cu and Fe ordering and crystal structure of NaFeo.56Cuo.44As, it is important to 
determine the magnetic states of Cu and Fe. Assuming that the covalent As-As bonding is not possible because of 
the large distance between the neighboring As ions [see Fig. 1(b)], As ions in NaFei-^Cu^As can only be in the As^“ 
state. From XAS and RIXS experiments on NaFeo.56Cuo.44As [Supplementry Note 2 and Supplementary Figure 3], 
we conclude that Cu in NaFei_xCua;As is in the nonmagnetic Cu^"*" configuration. As a consequence, Fe should be in 
the Fe^"*" 3c?® state since Na can only be in the Na^’*' state. 

Assuming Fe in NaFeo.56Cuo.44As is indeed in the Fe®+ 3d® state, we can determine the magnetic structure of 
the system by systematically measuring magnetic peaks at different wave vectors [Supplementary Note 3 and Sup¬ 
plementary Figures 4-5] and using neutron polarization analysis [Supplementary Note 3 and Supplementary Figure 
6]. By polarizing neutrons along the direction of the momentum transfer Q, neutron spin-flip (SF) scattering is 
sensitive to the magnetic components perpendicular to Q, whereas the non-spin-hip (NSF) scattering probes pure 
nuclear scattering 31, 3^. Figure 2(g) shows SF and NSF scans along the [IF, 0, 0.5]/[0, AT, 0.5]direction at T = 2 K 
and 240 K. Inspection of the data reveals clear magnetic scattering at 2 K that disappears at 240 K, in addition to 
the temperature independent NSF nuclear super-lattice rehection. Figure 2(h) shows the SF and NSF scans along 
the [id, id, 0] direction at 2 K and 240 K. While the SF scattering shows a clear peak at 2 K that disappears on 
warming to 240 K, the NSF scattering shows no evidence of the super-lattice peaks. From these results and from 
systematic determination of magnetic and super-lattice scattering intensity at different wave vectors, we conclude 
that NaFeo.56Cuo.44As forms a collinear magnetic structure with moments aligned along the a axis as shown in the 
right inset in Fig. 1(g). The ordered magnetic moment is 1.12 ±0.09 /is/Fe at 4 K determined from neutron powder 
diffraction [Supplementary Fig. 2(e)]. For this magnetic structure, magnetic peaks are expected and observed at 
(0,1,0.5) and (1,1,0) as shown in Fig. 2(g)-(h) [Supplementary Figure 7]. The magnetic peaks at (0,1,0.5) overlap 
with super-lattice peaks at (1,0,0.5) due to twinning. 


Doping dependence of magnetic order in NaFei_a:Cui As 

To determine the Cu-doping dependence of the NaFei_2;Cu2;As phase diagram, we carried out additional mea¬ 
surements on single crystals of NaFei_xCua;As with x = 0.18 and 0.39. Figures 3(a) and 3(b) compare the wave 
vector scans along the [IF, 0, 0.5]/[0, A", 0.5] and [1, 0, A]/[0,1, A] directions for x = 0.18,0.39, and 0.44. In each case, 
the scattering intensity is normalized to the (2,0,0) nuclear Bragg peak. With increasing Cu-doping, the scattering 
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FIG. 4: Summary of the schematic Cu and Fe electronic states and parameter regimes of a Mott insulator in 
the phase diagram. Schematic of the electronic states of a Cu ion alongside those of an Fe ion, shown in terms of the atomic 
levels (a) and the electronic density of states (b). The relative potential difference associated with the two ions is specified 
in terms of an energy shift, A, as illustrated in (b). (c) The ground-state phase diagram of a multiorbital Hubbard model 
for NaFe 0 .sCu 0 . 5 As. UosM and Umt refer to the critical U values for an orbital-selective Mott (OSM) transition and a Mott 
transition (MT), respectively into an orbital-selective Mott phase (OSMP) and a Mott insulator. The shaded region shows the 
physical parameter regime as determined by comparing the theoretically calculated bandwidth renormalization factors of the 
three t 2 g orbitals with those determined by the ARPES results [Supplementary Note 5]. 


profile becomes narrower and stronger, changing from a broad peak indicative of the short-range magnetic order 
in a; = 0.18,0.39 to an essentially instrument resolution limited peak with long-range magnetic order at a; = 0.44. 
Figure 3(c) shows the doping evolution of the spin-spin correlation length ^ in NaFei_a;Cua,As, suggesting that the 
increasing spin correlations in NaFei-ajCu^As are related to the increasing Cu-doping and the concomitant increase 
in correlation length of Fe and Cu ordering [Supplementary Note 3 and Supplementary Figure 5(f)]. The evolution 
of the temperature dependent magnetic order parameter with Cu doping is shown in Fig. 3(d). For x well below 0.5, 
the magnetic transition is gradual and spin-glass-like. However with x approaching 0.5, the transition at Tn becomes 
more well-defined. In the undoped state, NaFeAs has a small ordered moment of 0.17 ± 0.03 Mn/Fe [^ . Upon 
small Cu-doping, superconductivity is induced at a; = 0.02 and the static AF order is suppressed [^. With further 
Cu doping, the system becomes an AF insulator, where the ordered moment reaches ~1.1 ^g/Fe at a; = 0.44 (the 
ordered moment per Fe site was determined in the structure for NaFeo.sCuo.sAs [Fig. 1(d)]). Since the iron moment 
in NaFei-xCu^As increases with increasing Cu-doping, our determined moment is therefore a lower bound for the 
ordered moment for Fe in the ideal NaFeo.sCuo.sAs, which for Fe 3d® can be in either S = 3/2 or S' = 5/2 spin state. 


DISCUSSION 


The behavior in NaFei_a;Cua;As is entirely different from the bipartite magnetic parent phases seen in the iron 
oxypnictide superconductor LaFeAsOi-xHa,, where magnetic parent phases on both sides of the superconducting 
dome are metallic antiferromagnets (33| . In Cu-doped Fei-a^Cu^^Se, a metal-insulator transition has been observed 
around 4% Cu-doping, and a localized moment with spin glass behavior is found near x k, 0.12 [^[s^. Although the 
density functional theory (DFT) calculations suggest that Cu occurs in the 3d^® configuration and the metal-insulator 
transition is by a disorder induced Anderson localization in Fei_a;Cua;Se [s^, the 20-30% solubility limit of the system 
34| means it is unclear whether Fei_a;Cua;Se is also an AF insulator at x ~ 0.5. 

It is useful to compare our experimental findings with the results of the band structure calculations based on DFT 
[Methods section, Supplementary Note 4 and Supplementary Figures 8-10]. Both the DFT and DFT-bC/ calculations 
clearly predict the paramagnetic phase of NaFeo. 5 Cuo. 5 As to be a metal [Supplementary Figure 9]. By contrast, our 
measurements have shown that the insulating behavior of the resistivity persists above the Neel temperature [Fig. 
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1(f)], implying that NaFeo. 5 Cuo. 5 As is a Mott insulator. This is also consistent with STM measurements on lower 
Cu-doping NaFei_a;Cua;As where the overall line shape of the electronic spectrum at a; = 0.3 is similar to those of 
lightly electron-doped copper oxides close to the parent Mott insulator |21|. 

To understand the origin of the Mott insulating behavior, we address how Cu-doping affects the strength of electron 
correlations [Methods section. Supplementary Note 5 and Supplementary Figures 11-12]. Our starting point is that 
the local (ionization) potential difference between the Fe and Cu ions, as illustrated in Fig. 4(a) and described by the 
energy shift A in Fig. 4(b), will suppress the hopping integral between Fe and Cu sites. This causes a reduction in 
the kinetic energy or, equivalently, the electron bandwidth. The ratio of the on-site Coulomb repulsion U (including 
the Hund’s coupling Jh) relative to the bandwidth will effectively increase, even if the raw values of U and Jh remain 
the same as in pure NaFeAs. This would enhance the tendency of electron localization even without considering the 
disorder effects with increasing x 37|, l38 j . 

The case of x = 0.5 allows a detailed theoretical analysis. Here, the Cu and Fe ions in NaFei_a,Cua,As have real 
space ordering, as discussed above. We study the metal-insulator transition in multi-orbital Hubbard models for 
NaFeo. 5 Cuo. 5 As via the U{1) slave-spin mean-field theory [s^. We use the tight-binding parameters for NaFeAs, and 
consider the limit of a large local potential difference between Cu and Fe, which transfers charge from Fe to Cu to 
Cu^+ (n = 10). The fully occupied Cu 3d shell makes the Cu ions essentially as vacancies. Reduction of the Fe 3d 
electron bandwidth is caused by this kinetic blocking effect [Supplementary Figure 11]. The resulting ground state 
phase diagram is shown in Fig. 4(c). For realistic parameters, illustrated by the shaded region, a Mott localization 
takes place. Our understanding is in line with the general theoretical identification of a Mott insulating phase in 
an overall phase diagram [^ . which takes into account a kinetic-energy-reduction induced increase of the effective 
interactions and a decrease of the Sd-electron filling from 6 per Fe^^ to 5 per Fe^^ 

Our work uncovers a Mott insulator, NaFen. 5 Cun 5 AS, and provides its understanding in an overall phase diagram 
of both bandwidth and electron-filling controls [30|. Our results suggest that the electron correlations of the iron 
pnictides, while weaker than those of the iron chalcogenides, are sufficiently strong to place these materials in proximity 
to the Mott localization. This finding makes it natural that the electron correlations and the associated bad-metal 
behavior and magnetism induce a similarly high superconducting transition temperature in the iron pnictides as in 
the iron chalcogenide family. Such a commonality, in spite of the very different chemical composition and electronic 
structure of these two broad classes of materials, introduces considerable simplicity in the quest for a unified framework 
of the iron-based superconductivity. More generally, the proximity to the Mott transition links the superconductivity 
of the iron pnictides to that arising in the copper oxides, organic charge-transfer salts [s^ and alkali-doped fullerides 
0 , and suggests that the same framework may apply to all these strongly correlated electronic systems. 

Note added after submission: Very recently, angle-resolved pho toemission spectroscopy measurements on 
NaFeo. 56 Cuo. 44 As samples also confirmed its Mott insulating nature M|. 


METHODS 


Sample preparation and experimental details 

NaFei_a:Cua;As single crystals were grown by the self-flux method using the same growth procedure as for 
NaFei_a,Coa;As described in earlier work [^. The Cu doping levels used in this paper were determined by 
inductively coupled plasma (ICP) atomic-emission spectroscopy. Samples with nominal Cu concentrations of x = 
2%, 20%, 50%, 75% and 90% were prepared, resulting in actual Cu concentrations of x = 1.6%, 18.4(0.4)%, 
38.9(2.8)%, 44.2(1.7)% and 48.4(3.4)%. For each nominal doping except x = 2%, 5-6 samples were measured and 
the standard deviation in these measurements are taken as the uncertainty of the actual concentrations. For 
simplicity, the actual concentrations are noted as x = 1.6%, 18%, 39%, 44% and 48% in the rest of the paper. This 
suggests that the solubility limit of Cu in NaFei-xCu^As single crystals by our growth method is ~ 50%. 

For resistivity measurements, samples were mounted onto a resistivity puck inside an Ar filled glove box by the 
four-probe method and covered in Apiezon N grease. The prepared puck is then transferred in an Ar sealed bottle 
and is only briefly exposed in air while being loaded into a physical property measurement system (PPMS) for 
measurements. After measurements no visual deterioration of the samples were seen under a microscope. 







For single crystal elastic neutron scattering measurements, samples were covered with a hydrogen-free glue and then 
stored in a vacuum bottle at all times except during sample loading before the neutron scattering experiment. While 
our largest crystals grown are up to 2 grams, we typically used thin plate-like samples around 5x5 mm^ in size with 
mass ssO.2 grams for elastic neutron scattering measurements. Measurement of the sample with x = 0.18 was carried 
out on the N5 triple-axis spectrometer at the Canadian Neutron Beam Center (CNBC), Chalk River Laboratories. 
Pyrolitic graphite (PC) monochromator and analyzer {E^ = 14.56 meV) were used with none-36'-sample-33'-144' 
collimation setup. A PG filter was placed after the sample to eliminate contamination from higher order neutron 
wavelengths. The experiment on a; = 0.39 and 0.44 samples were done on the HB-3 triple-axis spectrometer and the 
HB-IA fixed-incident-energy triple axis spectrometer, respectively, at the High Flux Isotope Reactor, Oak Ridge 
National Laboratory. On HB-3, a PG monochromator {E^ = 14.7 meV), analyzer and two PG filters (one before and 
the other after the sample) were used with collimation setup 48'-60'-sample-80'-120'. HB-IA uses 2 PG 
monochromators {Ei = 14.6 meV) and 2 PG filters (mounted before and after second monochromator) resulting in 
negligible higher order wavelengths neutrons in the incident beam. A PG analyzer is placed after the sample. The 
collimation used was 48'-48'-sample-40'-68'. In all cases above, the samples were aligned in [iJ, 0, L] scattering plane. 
The polarized single crystal neutron scattering measurements were carried out on the G5 polarized-beam triple-axis 
spectrometer at GNBC, Ghalk River Laboratories. The neutron beams were polarized with Heusler (1,1,1) crystals 
with a vertically focusing monochromator and a flat analyzer (Af = 13.70 meV). A PG filter was placed after the 
sample and none-48'-51'-144' collimation was used. Neutron polarization was maintained by using permanent 
magnet guide fields. Mezei flippers were placed before the sample to allow the measurement of neutron SF and NSF 
scattering cross-sections. A 5-coil Helmholtz assembly was used to control the neutron spin-orientation at the 
sample position by producing a magnetic field of the order of 10 G. The orientation of the magnetic field at the 
sample position was automatically adjusted to allow the measurements to be performed for the neutron spin to be 
parallel or perpendicular to the momentum transfer. The flipping ratio, defined as the ratio of nuclear Bragg peak 
intensities in NSF and SF channels, was measured to be about 10:1 for various field configurations. The sample was 
studied in both [iJ, 0,L] and [iJ, iJ, L] scattering planes. 

The XAS and RIXS measurements on the x = 0.44 samples were carried out at the Advanced Resonant 
Spectroscopy (ADRESS) beam-line of the Swiss Light Source, Paul Scherrer Institut, Switzerland [i^. Samples were 
cleaved in-situ and measured in a vacuum better than 1 x 10“^*^ mbar. X-ray absorption was measured in total 
electron yield mode by recording the drain current from the samples and in total fluorescence yield with a 
photodiode. Linearly polarized X-rays with E = 931.8 eV resonant at the Cu edge were used for the RIXS 
measurements. The total momentum transfer was kept constant, but the component of momentum transfer in the 
ab plane has been varied. 

For neutron powder diffraction (NPD) measurements, the samples were ground from 2 grams {x = 0.016, 0.18, 0.39, 
0.44) of single crystals and sealed in vanadium sample cans inside a He filled glove box. NPD measurements were 
carried out at room temperature (300 K) on the BTl high resolution powder diffractometer at NIST Genter for 
Neutron Research. The Ge(3,1,1) monochromator was used to yield the highest neutron intensity and best 
resolution at low scattering angles. For the x = 0.44 sample, measurement at 4 K using a Gu(3,1,1) monochromator 
was also carried out. 

For the single crystal neutron diffraction measurement, a sample with x = 0.44 (25 mg) was used. The experiment 
was carried out at the four-circle diffractometer HB-3A at the High Flux Isotope Reactor, Oak Ridge National 
Laboratory. The data was measured at 250 K with neutron wavelength of 1.003 A from a bent Si(3, 3,1) 
monochromator using an Anger camera detector. Bragg peaks associated with the NaFeAs structure were measured 
with 1 second per point by carrying out rocking scans. Super-lattice peaks were measured with 10 minutes per point 
at each position. 


DFT-based electronic structure calculations 

In order to gain insight into the insulating nature of Cu-doped NaFeAs, we have performed a series of DFT based 
electronic structure calculations. The electronic band structure calculations have been performed using a full 
potential linear augmented plane wave method (LAPW) as implemented in the WIEN2K package [i^. The 
exchange-correlations function was taken within the generalized gradient approximation (GGA) in the 
parameterization of Perdew, Burke, and Ernzerhof (PBE) [4^. For the atomic spheres, the muffin-tin radii (Amt) 
were chosen to be 2.5 bohr for Na sites, 2.38 bohr for Fe and Cu, and 2.26 bohr for As sites, respectively. The 
number of plane waves was limited by a cutoff parameter (i?MT x ATmax) = 7.5. To account for the on-site 
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correlations in the form of the Hubbard parameter U, the GGA 
used 


U approach in a Hartree-Fock-like scheme was 


47l| with the value U = 3.15 eV as calculated in LiFeAs using the constrained-RPA approach |48|. 


While the neutron diffraction studies were performed on the a; = 0.44 sample, it is much easier theoretically to study 
the 50% Gu doped NaFeo.sGuo.sAs. To model the crystal structure, the experimentally determined lattice 
parameters and atomic coordinates of NaFeo, 56 Cuo. 44 As at 4 K were used [Supplementary Table 3]. The difference 
between the structure at a; = 0.44 and x = 0.5 was deemed inessential in view of the weak dependence of the lattice 
parameters on Cu concentration x [Supplementary Table 3]. While the experimental structure can be fitted with a 
symmetry group P4/nmm with disordered Gu/Fe, the high-resolution TEM measurements found that the 
super-lattice peaks indicate a stripe-like ordered pattern [Fig. Id] consistent with the space group Iham. We have 
therefore used the Ibam space group, with 4 formula units per unit cell, to perform the ab initio calculations. 
Interestingly, the results of the calculations turn out to be insensitive to the details of the atomic arrangement: for 
instance, arranging Gu/Fe atoms in the checkerboard pattern results in a very similar density of states (DOS) as 
that obtained for the stripe-like atomic configuration. We also note that the adopted ordered arrangement of Fe/Gu 
sites neglects the effects of disorder and Anderson localization, which may be important at low doping levels x < 0.5. 
Those effects are however beyond the grasp of the first principles DFT calculations. 


U{1) slave-spin mean-field theory 


To take into account the correlation effects beyond the Hatree-Fock level, we further study the metal-insulator 
transition in multi-orbital Hubbard models for both NaFeAs and NaFeo. 5 Guo. 5 As via the t/(l) slave-spin mean-field 
theory (49l [. The model Hamiltonian contains a tight-binding part and a local interaction part. Detailed description 
of the model is given in Ref. [d^. The tight-binding parameters for NaFeAs are taken from Ref. and the 
electron density is fixed to n = 6 per Fe. As for NaFeo. 5 Guo. 5 As, the very large local potential difference between Gu 
and Fe, as estimated from our and previous [37| DFT calculations, leads to charge transferring from Fe to Cu. As a 
result, Fe is effectively hole doped when the Cu doping concentration is close to 0.5. This picture is supported by 
our experimental results which suggest that the Cu is in a 3(i^*^ configuration with n = 10, and Fe configuration is 
3(i® with n = 5. As a zeroth-order approximation to the fully occupied Cu 3d shell, we treat Cu ions as vacancies. 
They organize themselves in a columnar-order pattern in NaFeo. 5 Cuo. 5 As, as shown in Fig. 1(d). 


Data availability 

The data that support the findings of this study are available from the corresponding authors upon request. 
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